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摘要 :【 目的 】 推 测 桔 小 实 晶 Bactrocera dorsalis 在 中 国 的 扩散 路 径 和 新 发 生地 区 的 入 侵 来 源 。[【 方 法 ] 本 研究 测序 
获得 来 自 中 国 ,泰国 日 本 ,老挝 .孟加拉 国 和 美国 等 地 31 个 种 群 的 192 头 桔 小 实 晶 个 体 的 COI 序列 (1 496 bp, 
COI 基因 全 长 97.3% ) ,并 以 软件 DnaSP 5.0, MEGA 6.0 和 Arlequin 3.5 等 完成 各 种 群 的 遗传 多 样 性 、 种 群 间 的 遗 
传 分 化 以 及 单 倍 型 分 析 。【 结 果 】 所 测 31 个 桔 小 实 晶 种 群 总 体 表 现 出 较 高 水 平 的 核 苷 酸 多 样 性 (0. 00663 ) 和 高 水 
平 的 单 倍 型 多 样 性 (0.98069) 。 以 Fr- 统计 法 度量 种 群 间 遗 传 分 化 程度 , 结果 显示 中 国 桔 小 实 晶 地 理 种群 间 遗传 分 
化 较 弱 ,中 国 种 群 与 泰国 .日 本 ,老挝 .和 孟加拉 国美 国 种 群 间 的 遗传 分 化 程度 不 同 , 中 国 种 群 与 美国 种 群 及 日 本 种 
群 的 遗传 分 化 最 大 。 而 Mantel 检验 发 现 ,中 国泰 国 日 本 ,老挝 .孟加拉 国 和 美国 种 群 间 的 遗传 分 化 与 空间 距离 有 
关 (R=0.670, P <0.0001), 中 国 种 群 间 的 遗传 分 化 不 是 地 理 隔离 所 造成 的 (R=0.038, P=0.534)。 中 国 种 群 与 
泰国 日 本 ,老挝 .孟加拉 国 和 美国 种 群 间 不 存在 共享 单 倍 型 。 根 据 种 群 系统 发 育 树 , 可 将 中 国 原 发 生地 区 种 群 划 
分 为 西南 .东南 两 大 分 支 。 中 性 检验 和 错 配 分 析 结 果 均 表明 桔 小 实 晶 曾 存在 大 规模 的 扩张 。【[ 结论] 村 小 实 蝇 以 东 
南 地 区 和 西南 地 区 为 源头 向 中 国内 陆 扩 散 ,其 中 广东 、 福 建 . 广 西 和 贵州 种 群 为 中 国内 陆 种 群 较为 有 影响 力 的 源 
头 。 根 据 遗 传 多 样 性 .遗传 分 化 . 单 倍 型 分 析 ,推测 新 发 生地 区 桔 小 实 晶 的 来 源 , 例如 安徽 合肥 桔 小 实 晶 种 群 主要 
来 源 于 福建 长 乐 广东 珠海 和 上 海 。 

关键 词 : RE СО] 基因 ; 种 群 遗传 分 化 ; 遗传 结构 ; 单 倍 型 ; 遗传 多 样 性 ; 种 群 扩张 
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Population genetic differentiation analysis of the oriental fruit fly, 
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Abstract: [ Aim] This study aims to infer the invasion routes of the oriental fruit fly, Bactrocera dorsalis 
Hendel in China and the invasion source of its recent introduction region. | Methods] mtDNA COI 
(1 496 bp, 97.396 of full length) of 192 individuals of 31 populations from China, Thailand, Japan, 
Laos, Bangladesh and USA was sequenced. The genetic diversity of each population and genetic 
differentiation between populations were analyzed by DnaSP 5.0, MEGA 6.0 and Arlequin 3. 5, as well 
as haloptype analysis. [ Results] High nucleotide diversity ( 0. 00663 ) and haplotype diversity 
(0.98069) were detected overall from the 31 populations of B. dorsalis. Genetic differentiation between 
populations analyzed using F-statistics method revealed that lower genetic differentiation were detected 
between populations of China, different levels of genetic differentiation were detected between populations 
from China and populations from Thailand, Japan, Laos, Bangladesh and USA, and the highest genetic 
differentiation was detected between populations from China and populations from USA and Japan. Mantel 
test showed that genetic differentiation among 31 populations from China, Thailand, Japan, Laos, 
Bangladesh and USA bore a relation to spatial distance (А 20.670, Р <0. 0001) , but spatial distance was 
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not responsible for genetic differentiation among populations from China (Ё 20.038, P 20.534). There 


was no shared haloptype among populations from China, Thailand, Japan, Laos, Bangladesh and USA. 


The phylogenetic tree showed that original populations from China clustered into two groups, i. e., the 


southwest group and the southeast group. Neutrality test and mismatch distribution analysis revealed that B. 


dorsalis populations had experienced a large-scale expansion. [Conclusion] It is inferred that there are two 


main invasion routes to inland China, one is from Southeast China and the other from Southwest China, and 


Guangdong, Fujian, Guangxi and Guizhou are the influential invasion sources. Based on the results of 


genetic diversity, genetic differentiation and haplotype network, it is inferred that the main sources of В. 


dorsalis introduced to Anhui ( Hefei) are from Fujian ( Changle) , Guangdong (Zhuhai) and Shanghai. 


Key words: Bactrocera dorsalis; COI gene; population genetic differentiation; genetic structure; 


haloptype; genetic diversity; population expansion 


吉 小 实则 Bactrocera dorsalis Hendel, XZ RIR 
SCH, Jm OH H (Diptera), Sz й F} ( Tephritidae ) 
(Hardy, 1969; 黄 德 聪 和 林 奇 力 ，1984 ) ， 为 害 柑橘 、 
Е У Е е ЕЕ 
KR, 是 热带 和 亚热带 水 果 的 主要 害虫 之 一 ( 李 红 旭 
等 , 2000; WEISE, 2008), 我 国 将 包括 桔 小 实 晶 在 内 
的 果实 蝇 属 Bactrocera spp. 均 列 为 进 境 植物 检疫 性 有 
害 生物 (中 华人 民 共 和 国 农业 部 , 2007). 















































江苏 .上海 浙江 等 新 发 生地 区 , 桔 小 实 蝇 的 发 生 面 
积 不 断 扩大 ,所 造成 的 经 济 损失 也 不 容 忽视 (吴广 
超 等 , 2007; 祁 力 言 等 , 2008; KERS, 2008; #0 
尧 平 和 赵 军 ,2011 ) 。 此 外 , 桔 小 实 晶 直 接 影响 水 
果 和 芯 菜 的 国际 贸易 , 许多 国家 对 我 国 桔 小 实 蝇 疫 
区 的 果蔬 有 严格 的 植物 检疫 准 入 要 求 ( 李 志 红 等 ， 
2002; WRR MSh Н ВЕ, 2005; THE, 2008 ) 。 
此 , 分 析 桂 小 实 蝇 种 群 间 的 遗传 进化 关系 ,推断 新 



































我 国 于 1911 年 在 台湾 首次 记录 桂 小 实 蝇 ( 李 文 
4€, 1988), 1934 年 在 海南 发 现 其 为 害 ( Wang, 
1996; RIPRESE, 2011), 20 20 50 年 代 后 , ЖУ 
实 蝇 在 广东 、 云 南 、 广 西 , 福 建 以 及 贵州 等 地 发 生 
(Wang, 1996; 林 振 基 等 , 1998; Б, 1998; ХІ 
建 宏 等 , 2010)。1980 年 , 四川 攀 枝 地 区 发 现 桂 小 
K OKE ZR, 1994) ,其 分 布 范围 持 
续 向 北 推进 ,陆续 在 江西 ,湖北 、 浙 江 、 上 海 、 江 苏 、 
安徽 和 重庆 等 地 区 发 现 其 为 害 ( 周 国 梁 等 ,2006; 
神力 言 等 , 2008; 赵 琳 等 , 2008; HEZ, 2010; 
万 宣 伍 等 ，2010; 王 永 模 等 ，2011; 能 桂 和 等 ， 
2011)， 目 前 该 虫 的 分 布 北 界 已 达 32?"2'N 的 无 锡 
( 祁 力 言 等 , 2008; 万 宣 伍 , 2012), 

随 着 桔 小 实 晶 分 布 范 围 的 扩大 和 发 生 量 的 增 
K, 桔 小 实 蝇 造 成 的 经 济 损失 不 断 提 高 。 在 厦门 ， 
吉 小 实 蝇 为 害 番 橄 榄 枇杷、 人 心 果 等 ,其 受害 率 在 
8096 以 上 , 厦门 岛 南 部 部 分 瓜 果 甚至 出 现 绝收 现象 
( 张 清 源 和 林 向 阳 ，1988; 林 振 基 等 ,1998 ) 。 在 云 
南 地 区 , 桔 小 实 晶 的 为 害 也 十 分 严重 , 水 果 的 受害 
率 曾 可 达 90% 以 上 ( 张 智 英 等 ,1995; 李 红 旭 等 ， 
2000) 2005 4E, 桔 小 实 晶 对 福建 省 造成 的 总 经 济 
损失 约 为 1.21 亿 元 ( 宁 昭 玉 , 2008). 2011 4E, ffi 
小 实 蝇 在 广东 省 的 发 生 面积 达到 240 khim, 包括 村 
小 实 蝇 在 内 的 3 种 实 蝇 对 广东 造成 的 经 济 损失 总 
估计 在 33.67 ~ 129.87 亿 元 ( 马 兴 莉 等 , 2013) 。 在 




























































































发 生地 区 村 小 实 蝇 的 种 群 来 源 , 对 该 虫 的 防 控 工作 
具有 重要 意义 。 

线粒体 基因 因 其 母系 遗传 .几乎 不 发 生 重 组 ЛЕ 
细胞 中 拷贝 数量 高 等 特点 , 成 为 种 群 间 遗传 分 化 研 
究 的 主要 分 子 标记 ( 何 恒 果 ,2008; Rollins et al., 
2011)。 云 南 地 区 是 目前 村 小 实 蝇 种 群 间 遗传 分 化 
研究 最 多 的 地 区 ( 施 伟 和 叶 辉 ,2004; Shi et al., 
2005; Liu et al., 2007; 施 伟 和 叶 辉 ，2007; Shi et 
al., 2010) , 通过 对 一 个 地 区 的 村 小 实 蝇 种 群 间 的 
遗传 分 化 研究 , 可 判断 该 地 区 桔 小 实 晶 的 人 侵 来 
源 , 为 检疫 和 防 挖 提供 理论 支持 (万 宣 伍 等 ,2010; 
吴广 超 等 , 2012) 。 当 桔 小 实 蝎 种 群 间 遗 传 分 化 研 
究 扩大 至 全 球 范 围 时 ， 可 对 桔 小 实 蝇 的 起 源 .扩散 
路 径 进行 推测 (Aketarawong et al., 2007; FESE, 
2007; Wan et al., 2011, 2012; Li et al., 2012; Shi et 
al., 2012) , Li 4: (2012) 以 COI 序列 对 15 个 村 小 
实 蝇 种 群 进行 遗 传 分 化 分 析 , 结果 显示 桔 小 实 晶 起 
源 于 东南 亚 。Shi 等 (2012 ) 对 29 个 桔 小 实 晶 地 理 
种 进行 遗传 分 化 研究 , 结果 显示 村 小 实 蝇 可 能 起 源 
于 东南 亚 的 热带 地 区 和 中 国 的 南部 沿海 地 区 。 而 
Wan 等 (2011, 2012) 的 研究 结果 显示 , 桔 小 实 蝎 可 
能 起 源 于 中 国 台 湾 , 并 通过 东南 亚 地 区 和 中 国 东 南 
沿海 地 区 向 中 国 大 陆 扩散 。 然 而 , 目前 中 国 村 小 实 
蝇 地 理 种 群 遗传 分 化 研究 的 取样 范围 较 小 ， 且 取样 
地 点 较为 集中 ( 李 伟 丰 等 , 2007; Wan et al., 2011, 
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2012; Li et al., 2012; Shi et al., 2012), ВЖ, 本 研 
究 进一步 扩大 取样 范围 , 测评 全 国 现 有 分 布 记录 的 
省 级 行政 区 的 桔 小 实 晶 种 群 遗 传 分 化 情况 ,较为 全 
面 地 分 析 村 小 实 蝇 在 中 国 的 扩散 路 径 和 新 发 生地 区 
的 入 侵 来 源 。 

本 研究 系统 收集 了 国内 外 村 小 实 蝇 31 个 地 理 
种 群 共 192 头 成 虫 , 通过 DNA 测序 获得 线粒体 COI 
基因 ,对 我 国 不 同 地 理 种 群 的 桔 小 实 晶 基因 信息 、 
遗传 多 样 性 和 种 群 亲缘 关系 进行 了 研究 , 为 防 控 我 
国 材 小 实 蝇 的 持续 扩散 提供 依据 。 



































1 材料 与 方法 


1.1 і 

供 试 桔 小 实 蝇 共有 31 个 种 群 192 头 成 虫 标本 ， 
样品 来 源 包括 我 国 现 有 分 布 记录 的 省 级 行政 区 、 东 
南亚 和 美国 夏威夷 地 区 ( 表 1), 其 中 中 国 23 个 种 
WE, 泰国 4 个 种 群 ,日 本 老挝、 孟加拉 国 和 美国 各 
1 个 种 群 。 本 研究 根据 历史 记录 , 将 我 国 的 桔 小 实 
蝇 种 群 初步 划分 为 2 Ж. 原 发 生地 区 ( 台北、 香港 、 
澳门 .万 宁 、 深 圳 \ 珠 海 、 河 口 瑞丽、 动 腊 、 东兴 、\ 长 
乐 . 贯 阳 和 成 都 ) 种群 和 新 截获 或 发 生地 区 (合肥 、 
ЖЖ E BH ЕКЕЖ, Ж LE ZR ЛМ ЕНУ ВР FIN 
隆 ) 种 群 。 所 有 实 晶 成 虫 标本 均 由 深圳 出 人 境 检验 
检疫 局 动 植物 检验 检疫 技术 中 心 收集 和 鉴定 确认 ， 
均 用 无 水 酒精 浸泡 , 于 -20Y 保存 。 

1.2 DNA 提取 

桔 小 实 晶 单 头 个 体 总 ОМА 提取 以 及 质量 检测 方 
法 参照 Yu 等 (2005) 。DNA 样品 保存 于 -20% 备 用 。 
1.3 PCR 扩 增 及 序列 测定 

线粒体 СОГ 基因 (mtDNA COI) 由 UEA7/ 
UEAIO ( Lunt et al., 1996) 和 COIF/COIR ( Tan, 
2004)2 对 引物 进行 PCR 扩 增 , 最 后 经 拼接 获得 
СОТ 序列 1 496 bp, 引物 由 华 大 基因 合成 。 

PCR 反应 体系 为 25 pL, 包括 20 ~ 50 ng 模板 
DNA, 1 x PCR buffer (不 含 Mg ), 2. 5 mmol/L 
Mg'* , 0.2 mmol/L dNTPs, 1 U Taq 酶 ,上 下 游 引物 
40.2 mmol/L， 加 去 离子 水 至 总 体积 为 25 uL, 2 
对 引物 的 PCR. 扩 增 反应 条 件 为 ,引物 UEAT/ 
СЕА10; 94^C 预 变 性 4 шіп, 然后 运行 30 个 循环 ( 
个 循环 包括 950C 变性 40 s, 48% 退火 30 s, 72%C 延 
伸 1 min), 72°С 延伸 7 min。 引 物 COIR/COIF: 
94°С 预 变性 4 min, 然后 运行 35 个 循环 (每 个 循环 
包括 94°С 变性 40 s, 55% 退火 1 min, 72°С 延伸 2 















































шіп), 72°С 延伸 15 min, PCR 产物 于 1 % 琼脂 糖 凝 
胶 ( 加 EB) 上 检测 , 检测 结果 为 阳性 后 委托 华 大 基 
因 进 行 测序 , 测序 工作 在 ABI 3730 测序 仪 上 进行 ， 
所 有 序列 均 采 用 双向 测序 以 保证 序列 的 准确 性 。 
1.4 数据 分 析 

将 两 对 引物 扩 增 测序 所 得 序列 由 Seqman 
( DNAStar 软件 包 ) 进行 剪 切 和 拼接 ,人工 重复 校对 
3 次 , 在 软件 MEGA 6.0 中 (Tamura et al., 2013) 3t 
行 序列 比 对 。 多 态 性 位 点 数 、 简 约 信息 位 点 数 和 自 
裔 位 点 数 由 МЕСА 6.0 完成 。 桔 小 实 量 各 种 群 内 的 
变异 位 点 数 、 核 音 酸 多 样 性 ,平均 核 苷 酸 差异 . 单 倍 
型 数量 . 单 倍 型 多 样 性 在 DnaSP 5.0 软件 (Librado 
and Rozas, 2009) 中 计算 。 种 群 间 遗 传 分 化 系数 Fs 
由 软件 Arlequin 3. 5 ( Excoffier and Lischer, 2010) 计 
算 。 以 Excel 插件 XLSTAT( http ;//www. xlstat. com/ 
en/) 计算 两 方面 内 容 : a) 遗传 多 样 性 指数 与 样本 数 
量 之 间 的 Pearson 相关 系数 ; b) 完成 Mantel 检验 ， 
分 析 地 理 距 离 和 遗传 分 化 系数 间 的 相关 性 。 以 软件 
Migrate 3.6.4(Beerli апа Felsenstein, 2001; Beerli, 
2006) 计算 种 群 间 迁 移 率 ,其 参数 设 定 为 : 以 贝 叶 
斯 推理 (BI) 搜索 策略 运行 1 条 长 链 , 重复 搜索 
10 000 000 代 3 x, 27 10 000 代 老 化 样本 
(burn-in) 。 第 一 次 运行 的 9 及 M 的 初始 值 由 Fs 
估计 所 得 , 以 前 一 次 运行 的 结果 作为 后 一 次 的 初始 
值 , 重复 运行 Migrate 直至 数据 稳定 。 基 于 碱 基 转 
АЛП 6 Н] Kimura 双 参 数 模 型 ， 由 МЕСА 6.0 
(Tamura et al., 2013 ) 计 算 各 种 群 两 两 之 间 的 遗传 
ІН, 由 软件 Phylip 3. 69 ( Felsenstein, 2005) 基于 
遗传 距离 矩阵 ， 以 邻接 法 (NJ) 构建 种 群 系统 发 育 
Ж. ІН DnaSP 5.0 软件 (Librado and Rozas, 2009) 
进行 中 性 检验 和 错 配 分 布 检验 。 以 DnaSP 5.0 软件 
( Librado and Rozas, 2009) 对 COI 序列 进行 单 倍 型 
判定 , 即 碱 基 排 列 完全 相同 的 序列 为 同一 单 倍 型 。 
由 Network 软件 (Bandelt ег al., 1999; Polzin and 
Daneshmand , 2003) 构建 单 倍 型 网 络 图 ， 对 图 进行 
调整 和 校对 。 清 到 和 冲绳 种 群 因 样 本 数量 为 1, 无 
法 计算 遗传 多 样 性 指数 。 

本 研究 应 用 1 496 bp СОТ 序列 研究 村 小 实 蝇 种 
群 的 遗传 分 化 。 目 前 , 应 用 于 村 小 实 蝇 种 群 遗传 分 
化 研究 的 СОТ 序列 长 度 主 要 为 500 ~ 800 Бр( ЖИЛ, 
2006; Wan et al., 2011; Shi et al., 2012) ， 为 评价 
1 496 bp COI 序列 在 桔 小 实 蝇 种 群 遗 传 分 化 研究 中 
的 效用 , 根据 Wan 2 (2012) 研究 中 COL 序列 
(JN643923) , 将 本 研究 所 获得 1 496 bp COL 序列 进 
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ЖІ 桔 小 实 蝇 样 本 信息 和 线粒体 CO 基因 序列 号 


Table 1 Sample information of Bactrocera dorsalis and GenBank accession number of mitochondrial COI genes 

































































种 群 代码 采样 地 "m 经 度 样本 数 采样 日 期 GenBank 登录 号 
Population code Collecting locality Latitude Longitude Number id Collecuing dates GenBank acvession 
specimens (year. month) number 
DCHF 中 国安 徽 合 肥 Hefei, Anhui, China 117?13'E 31?49'N 5 2012.10 KF998619 — KF998623 
DCJS 中 国 江 苏 苏州 Suzhou, Jiangsu, China 120°35'Е 31?17'N 22 2012.10 KF998657 — KF998678 
DCSH 中 国 上 海 Shanghai, China 121°24'Е 31?13'N 6 2012.9 KF998702 — KF998707 
DCZJ 中 国 浙江 富阳 Fuyang, Zhejiang, China 119 °57'Е 30?03'N 10 2012.10 KF998777 - KF998786 
DCSC 中 国 四 川 成 都 Chengdu, Sichuan, China 104°03'Е 30?34'N 5 2012.9 KF998697 — KF998701 
DCWH 中 国 湖北 武汉 Wuhan, Hubei, China 114°17'Е 30?35'N 5 2012.10 KF998740 — KF998744 
DCCQ 中 国 重 庆 武 隆 Wulong, Chongqing, China 106°32'Е 29?33'N 5 2012.10 KF998597 — KF998601 
DCZM 中 国 西藏 樟 木 Zhangmu, Tibet, China" 859557Е 27?50'N 10 2013.9 KF998749 — KF998758 
DCYD 中 国 西藏 亚 东 Yadong, Tibet, China * 88°53Е 279287 4 2012.10 КЕ998745 - КЕ998748 
DCHY 中 国 湖南 衡阳 Hengyang, Hunan, China 112°33'Е 26?53'N 16 2012.8 KF998640 - KF998655 
DCGZ 中 国 贵 州 贵阳 Guiyang, Guizhou, China 106°37'Е 26?38'N 5 2011.10 KF998614 — KF998618 
DCFJ 中 国 福建 长 乐 Changle, Fujian, China 119°30'Е 25957М 4 2012.8 КЕ998610 - КЕ998613 
DCJX 中 国 江 西 竟 州 Ganzhou, Jiangxi, China 114°55'Е 25?49'N 5 2011.8 KF998679 - КЕ998683 
DCTW 中 国 台湾 台北 Taipei, Taiwan, China 121°33'Е 25?05'N 5 2012.8 КЕ998735 — КЕ998739 
DCYH 国 云南 河口 Некои, Yunnan, China 103°55'Е 22?3]'N 4 2012.10 KF998759 — KF998762 
DCYM 中 国 云南 动 腊 Mengla, Yunnan, China — 101?33'E 21?27'N 5 2012.10 KF998763 — KF998767 
DCYR 中 国 云南 瑞丽 Ruili, Yunnan, China 97°50'Е 24%07М 5 2012.10 KF998768 - KF998772 
DCSZ 中 国 广东 深圳 Shenzhen, Guangdong, China 114°02'Е 22?32'N 20 2012.9 KF998708 — KF998727 
DCZH 中 国 广东 珠海 Zhuhai, Guangdong, China 113°34'Е 22?16'N 4 2011.11 КЕ998773 — КЕ998776 
DCDX 中 国 广西 东兴 Dongxing, Guangxi, China 107°57'Е 21?32'N 8 2012.8 KF998602 — KF998609 
DCHK 中 国 香港 Hong Kong, China 114?06'E 22?23'N 5 2012.9 КЕ998624 - КЕ998628 
DCMC 中 国 澳 门 Macau，China 113°32'Е 22°11'\ 9 2012.9 KF998684 - КЕ998692 
DCHN 中 国 海南 万 宁 Wanning, Hainan, China 110°22'Е 18?47'N 5 2013.8 KF998629 — KF998633 
DLLN 老挝 琅 南 塔 Luang Namtha, Laos ** 101?09'E 20?55'N 4 2012.10 KF998693 - КЕ998696 
DTCM [8895 Chiang Mai, Thailand ** 98?59'E 18?47'N 1 2003.7 КЕ998728 
DTSI 泰国 四 色 菊 府 Sisaket, Thailand ** 104°18'Е 15?07'N 2 2013.6 KF998733 — KF998734 
DTPT 泰国 巴 知 他 尼 府 Pathum Thani, Thailand ** 100°31'Е 14°01' 2 2013.6 КЕ998729 – КЕ998730 
DTSA 泰国 沙 缴 府 Sa Kaew, Thailand ** 102*03'E 13?49'N 2 2013.6 KF998731 - KF998732 
DBDI 2. Е 88°37'Е 25?37'N 5 2013.9 KF998589 - KF998593 
DJNA 日 本 冲 蝇 Okinawa, Japan ** 127°40'Е 26?12'N 1 2004.7 КЕ998656 
DUHA 美国 夏威夷 Hawaii, USA ** 157952/7/ 219187 6 2004.7 KF998634 - КЕ998639 














* 口岸 检疫 截获 标本 Specimens from port quarantine interception; ** 国内 外 同行 赠送 Donated by domestic and foreign counterparts. 其 余 样 品 为 实 


晶 监 测 捕 获 标本 The rest of samples were obtained by monitoring. 


行 剪 切 ， 获 取 与 之 相应 的 短片 段 COI 序列 (561 
bp) ,计算 短片 段 СОГ 序列 的 多 态 性 位 点 数 .简约 信 
息 位 点 数 和 自 裔 位 点 数 等 序列 基本 数据 ,并 计算 单 
倍 型 数 、 核 苷 酸 多 样 性 .平均 核 苷 酸 差异 数 和 单 倍 型 
多 样 性 共 4 项 遗传 多 样 性 指数 。 








2 结果 


2.1 ҖЕ ЕНЕ COI 基因 的 РСЕ 扩 增 及 序 
列 变异 分 析 
通过 两 对 引物 UEAT/UEAIO 和 COIRZCOIF 对 村 








小 实 蝇 线 粒 体 COL 基因 进行 扩 增 , 经 拼接 校对 , 最 后 
获得 1 496 bp COI 序列 , 占 COI 全 长 的 97.3% 192 
条 СОГ 序列 已 全 部 提交 到 GenBank, 序列 号 见 表 1。 

对 192 Ж COI 基因 (1 496 bp ) 序 列 进 行 分 析 ， 
发 现 自 裔 位 点 83 A, HAE DNA 序列 的 5. 55% ; 
简约 信息 位 点 110 个 , 占 该 段 DNA 序列 的 7.35% ; 
总 多 态 性 位 点 193 个 ， 占 该 段 DNA 序列 的 
12.90% 。 其 中 含 2 个 突变 位 点 的 自 裔 位 点 80 个 ， 
简约 信息 位 点 105 个 ; 83 个 突变 位 点 的 自 裔 位 点 
3 个 , 简约 信息 位 点 5 个 , 自 裔 位 点 和 简约 信息 位 
点 均 未 发 现 4 个 突变 位 点 。 
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2.2 hKIR EE ЖЕНЕ 

Ш ЛУЗЕШ ЕЕЕ ЛЕ 3H АС, РЕВ ЖЖ 
酸 多 样 性 .平均 核 昔 酸 差异 数 和 单 倍 型 多 样 性 等 5 
项 遗传 多 样 性 指数 见 表 2。 村 小 实 蝇 31 个 种 群 总 
体 表 现 出 较 高 水 平 的 核 昔 酸 多 样 性 (0. 00663) 和 高 
水 平 的 单 倍 型 多 样 性 (0. 98069) , 各 种 群 的 遗传 多 
单 性 指数 表现 出 一 定 的 差异 。 另 外 , 各 种 群 的 样本 
数量 与 核 音 酸 多 样 性 平均 核 苷 酸 差异 数 . 单 倍 型 多 
样 性 的 Pearson 系数 分 别 为 -0.0578( 忆 =0.7663 ) , 
-0.0579( 忆 =0.7654) 和 - 0. 0194 (P =0. 9205), 
表明 种 群 间 的 遗传 多 样 性 差异 与 样本 数量 无 关 。 各 
种 群 的 核 背 酸 多 样 性 在 0.00089 -0.01738, 平均 核 
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苷 酸 差异 数 在 1.33333 ~ 26. 00000 之 间 , 其 中 美国 
夏威夷 种 群 核 背 酸 多 样 性 和 平均 核 音 酸 差异 数 最 
(Қ, 泰国 沙 缴 府 种群 核 并 酸 多 样 性 和 平均 核 昔 酸 差 
异 数 最 高 ; 单 倍 型 多 样 性 范围 是 0. 33333 ~ 
1.00000, 夏威夷 单 倍 型 多 样 性 最 低 ， 为 0. 33333, 

其 余 种 群 单 倍 型 多 样 性 均 较 高 , 其 中 62.19 的 种 群 
单 倍 型 多 样 性 为 1 00000, 4722 可 知 , 夏威夷 种 
群 的 遗传 多 样 性 较 低 , 泰国 种 群 遗 传 多 样 性 较 高 。 
中 国 部 分 新 截获 或 发 生地 区 种 群 遗传 多 样 性 高 于 原 
发 生地 区 种 群 的 遗传 多 样 性 ， 如 衡阳 种 群 和 亚 东 种 
群 的 核 苷 酸 多 样 性 和 平均 核 苷 酸 差异 数 大 于 除 东 兴 
种 群 外 的 原 发 生地 区 种 和 
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表 2， 桔 小 实 晶 各 种 群 遗 传 多 样 性 指数 和 单 倍 型 分 布 
Table 2 Genetic diversity index and haloptype distribution of Bactrocera dorsalis populations 
种 群 yo s m F га 单 倍 型 分 布 (频率 * ) 
Population Haloptype distribution ( frequency * ) 

DCCQ Э 17 0.00508 7.60000 1.00000 Н4, Н5, Нб, Н7, Н8 

рсрх 7 42 0.00869 13.00000 0. 96429 H9, H10, H11, H12, НІЗ, Н14(2), H15 

DCFJ 3 13 0.00479 7.16667 0.83333 H14(2), H16, H17 

DCGZ 5 16 0.00455 6.80000 1.00000 Н16, НІ8, НІ9, Н20, Н21 

DCHF 5 23 0.00668 10.00000 1.00000 H14, H22, H23, H24, H25 

DCHK 3 17 0.00535 8.00000 0.70000 НІ4(3), H26, H27 

DCHN 5 17 0.00508 7.60000 1.00000 H28, H29, H30, H31, H32 

DCHY 13 56 0.00730 10.91667 0. 96667 НАШ) Ну B33, Has, H33, M36, H37, Н38(3),„ 
H39, H40, H41, H42, H43 

DCJS 14 50 0.00624 9.33333 0.93506 HECA dz тана Нау Eos HAG: 

H47, H48, H49, H50, H51, H52, H53 

DCJX 3 11 0.00321 4.80000 0. 70000 H14, H31(3) , H54 

DCMC 9 36 0.00661 9.88889 1. 00000 H31, H55, Н56, Н57, H58, Н59, Н60, H61, H62 

DCSC 5 23 0.00682 10.20000 1.00000 H63, H64, H65, H66, H67 

DCSH 5 20 0.00579 8.66667 0. 93333 H14, H24, H68, Н69, H70(2) 

DCSZ 14 54 0.00660 9.86842 0. 93684 HIES y HAIZ) Аы H73, HTA, 

H75(2), H76, H77, H78, H79, H80, H81, H82 

DCTW 4 22 0.00628 9.40000 0.90000 H17, H31(2) , H83, H84 

DCWH 5 17 0.00521 7.80000 1.00000 Н14, H31, H85, H86, H87 

DCYD 3 16 0.00713 10.66667 1.00000 H88, H89, H90 

DCYH 4 20 0.00680 10.16667 1.00000 Н91, H92, H93, H94 

DCYM 5 22 0.00615 9.20000 1.00000 H95, H96, H97, H98, H99 

DCYR 5 20 0.00602 9.00000 1.00000 H100, H101, H102, H103, H104 

DCZH 4 14 0.00468 7.00000 1.00000 Н72, Н105, Н106, НІ07 

DCZJ 9 29 0.00532 7.95556 0.97778 H32(2) , H55, H72, H108, H109, H110, H111, H112, H113 

DCZM 10 24 0.00429 6.42222 1.00000 H7, H19, H104, H114, H115, H116, H117, H118, H119, H120 

DBDI 3 14 0.00624 9.33333 1. 00000 H1, H2, H3 

DJNA 1 N/A N/A N/A N/A Н121 

DLLN 4 19 0.00657 9.83333 1. 00000 Н122, Н123, Н124, Н125 

DTCM 1 N/A N/A N/A N/A H126 

DTPT 2 6 0.00401 6.00000 1.00000 Н127, Н128 

DTSA 2 26 0.01738 26.00000 1.00000 Н129, НІЗ0 

DTSI 2 12 0.00802 12.00000 1.00000 НІЗІ, H132 

DUHA 2 4 0.00089 1.33333 0. 33333 Н121(5), H133 

лк Е 


S. 变异 位 点 数 Number of variable sites; Н: 
数 Average number of nucleotide differences; На: FAAM LH 














单 倍 型 数量 Number of unique haplotypes; Pi: РЕ Nucleotide diversity; К; 
性 Haplotype diversity ; 


平均 核 并 酸 差异 
“ 单 倍 型 频率 默认 为 1, 不 具体 标明 Default value of haloptype 








frequency is 1 which is not displayed in detail; N/A; 单一 样本 , 不 适宜 计算 遗传 多 样 性 One sample only, genetic diversity is not applicable. 
































12 } E RF: 基于 线粒体 CO 基因 的 桔 小 实 晶 种 群 遗 传 分 化 研究 1429 














2.3 ”村 小 实 蝇 种 群 的 遗传 分 化 分 析 

采用 固定 指数 法 (或 称 Ff- 统计 ) 上 度量 种 群 间 的 
遗传 分 化 程度 ( Chakraborty and Danker-hopfe, 
1991), 240.00<„<0.05 时 , 遗传 分 化 程度 很 
Яя; 0.05 < Fs, «0. 15 时 ,遗传 分 化 程度 为 中 等 ; 
0. І5<Ғ, x0. 25 时 ， 遗传 分 化 程度 较 大 ; Ға > 
0.25 时 ,遗传 分 化 程度 很 大 (Rousset，1997; 万 宣 
(н, 2012) 31 个 村 小 实 蝇 种 群 间 的 遗传 分 化 系数 
见 表 3。 中 国 桔 小 实 蝇 种 群 间 遗传 分 化 程度 不 同 ， 
其 中 存在 基因 交流 的 占 42. 69% (Fo < 
-0.00018) , 遗传 分 化 程度 很 弱 、. 中 等 和 较 大 的 分 
别 占 41.11% ‚ 15. 42% 和 0. 79% (0. 00000 < Р... 
0.18632) ; 孟加拉 国 种 群 与 15 个 中 国 种 群 存在 基 
ЗЕЙ Рэг -0.00804)， 分 别 与 7 和 1 个 中 国 种 
群 存在 很 弱 .中 等 程度 的 遗传 分 化 (0. 00231 К< 
0.10654); 老挝 种 群 与 7 个 中 国 种 群 存在 基因 交流 
(Ес -0.00498) , 分 别 与 9,6 和 1 个 中 国 种 群 存 
在 很 弱 .中 等 和 较 大 程度 的 遗传 分 化 (0.01122 三 Fw 
«0. 20347) ; 泰国 种 群 与 中 国 种 群 间 的 遗传 分 化 程 
度 不 等 ,其 中 种 群 间 存 在 基因 交流 的 占 11. 96% 
(Fsrs - 0. 01215), 种 群 间 的 遗传 分 化 程度 为 很 
弱 、 中 等 、 较 大 和 最 大 的 分 别 占 11. 96% , 36. 96% , 
27.1796 和 11.96% (0. 0066-<Ғ;,<<0.45455); HÆ 
冲 强 种 群 与 中 国 种 群 间 存在 很 大 程度 遗传 分 化 






























































本 冲绳 、 美 国 夏威夷 种 群 间 存在 不 对 称 迁 移 ， 中 国 
原 发 生地 区 与 新 发 生地 区 种 群 间 存 在 不 对 称 迁移 。 
2.5 桔 小 实 蝇 种 群 的 单 倍 型 分 析 

192 条 桂 小 实 蝇 COI 序列 共 鉴 定 出 单 倍 型 133 
个 , 各 种 群 单 倍 型 的 组 成 见 表 2。 其 中 共享 单 倍 型 
16 个 , XR 117 个 (87. 97% ) 为 单个 种 群 所 特有 。 
除 美国 夏威夷 种 群 外 ， 独 享 单 倍 型 所 占 的 比率 以 及 
分 布 频 率 均 较 高 , 表明 这 些 地 理 种 群 之 间 已 呈现 出 
一 定 程度 的 遗传 分 化 。 中 国 种 群 与 东南 亚 、 美 国 夏 
威 夷 种 群 间 不 存在 共享 单 倍 型 。 

单 倍 型 网 络 图 (图 1) 总 体 呈 现 星 状 分 布 图 , 说 
明 桂 小 实 蝇 经 历 过 种 群 扩张 过 程 。 低 频率 单 倍 型 分 
布 在 单 倍 型 网 络 的 外 围 , 经 过 一 步 或 几 步 突变 与 高 
频率 单 倍 型 相连 。 与 中 国 原 发 生地 区 种 群 相 比 , 中 
国 新 发 生地 区 种 群 位 于 单 倍 型 网 络 的 外 围 。 单 倍 型 
Н1,Н2 和 НЗ 为 孟加拉 国 朗 布 尔 专区 种 群 特有 单 倍 
型 , 其 中 НЗ 经 两 步 突 变 由 H7( 武 隆 和 榜 木 ) 衍 生 ; 
泰国 巴 大 他 尼 府 种 群 特有 单 倍 型 H128 衍生 于 单 倍 
型 H19( 贵 阳 和 樟 木 ) ， 除 单 倍 型 H128 外 , RE A 
挝 ,日 本 种 群 与 中 国 种 群 无 衍生 关系 。 
26 ， 桔 小 实 蝇 种 群 系统 发 育 分 析 

种 群 系统 发 育 树 (图 2) 显示 : 孟加拉 国 朗 布尔 
专区 种 群 自 成 1 文 ,其 他 种 群 可 分 为 2 Ж. ЖЖ, 
澳门 .香港 .万 宁 、 成 都 . 动 脂 、 瑞 丽 、 老 挝 琅 南 塔 、 泰 

























































































(0.28767 < F,, «0. 69231); 美国 夏威夷 种 群 与 中 
国 种 群 间 存 在 很 大 程度 遗传 分 化 (0. 54899 < F or < 
0.81815) 。 总 体 上 ,中 国 种 群 与 孟加拉 国 、 老 挝 种 
群 的 遗传 分 化 程度 很 弱 , 与 泰国 种 群 间 遗 传 分 化 程 
度 中 等 或 较 大 , 与 日 本 、 美 国 种 群 间 遗 传 分 化 程度 
很 大 。 中 国 种 群 间 存在 基因 交流 或 遗传 分 化 程度 很 
Яя, 其 中 中 国 原 发 生地 区 中 , RE RKN KR 
阳 种 群 与 较 多 的 新 发 生地 区 种 群 存 在 基因 交流 。 

对 种 群 Rxx 值 和 空间 距离 两 个 矩阵 间 的 Mantel 
检验 , 结果 显示 : 全 部 种 群 的 两 个 矩阵 间 存 在 显著 
的 相关 性 (R=0.670, P <0.0001), 中 国 种 群 的 两 
个 和 矩阵 间 不 存在 显著 的 相关 性 (R=0.038, Р = 
0.534)。 表 明 中 国 、 泰 国 \ 日 本 、 老 挝 、 孟 加 拉 国 和 
夏威夷 种 群 间 的 遗传 分 化 与 空间 距离 有 关 , 中 国 种 
群 间 的 遗传 分 化 不 是 地 理 隔离 所 造成 的 。 
2.4 ”村 小 实 昌 种群 间 的 有 效 迁 移 率 

当 两 种 群 间 相互 的 有 效 迁 移 率 的 95% 置信 限 
不 重合 时 ， 即 认为 这 两 个 种 群 间 存在 不 对 称 迁移 
(万 宣 伍 , 2012)。 在 31 个 种 群 间 检 测 到 32 对 不 对 
称 迁 移 率 , 结果 如 表 4 所 示 。 中 国 与 泰国 老挝 日 










































































国 巴 吞 他 尼 府 以 及 泰国 四 色 菊 府 等 聚 为 分 支 1 
(clade 1) ,长乐 深圳、 珠海 . 闵 州 、 台 北 和 河口 等 聚 
在 分 支 2(clade 2), 新近 发 生地 区 种 群 在 2 个 分 支 
中 均 有 人 分布。 老挝、 美国 .泰国 和 日 本 种 群 均 位 于 分 
Xl, 

根据 种 群 系统 发 育 树 ,可 将 中 国 原 发 生地 区 种 
群 划分 为 两 大 分 支 : (1) 西南 分 支 : 东兴 、 贵 阳 、 澳 
门 . 香 港 、 万 宁 、 成 都 . 动 腊 和 瑞丽 种 群 ; (2) 东南 分 
ж: 台北 深圳、 珠海 .长 乐 和 河口 种 群 。 
2.7 桔 小 实 蝇 种 群 的 错 配 分 析 

31 个 桔 小 实 蝇 种 群 作 为 1 个 整体 进行 中 性 检 
验 和 错 配 分 布 分 析 ， 结 果 显 示 Тарта” s D 为 
-2.26055, HE МЕТЫ (Р <0.01); Fu's Fs 
JX - 169.055, 达到 显著 性 水 平 (P<0.02)。 错 配 
分 析 结 果 ( 3) 显示 为 单 峰 结构 。 中 性 检验 和 错 配 
分 析 结 果 均 证 明 桔 小 实 蝇 曾 存在 大 规模 的 扩张 。 

为 探究 村 小 实 晶 在 中 国 的 种 群 动态 , 以 23 个 
中 国 桔 小 实 晶 种 群 为 分 析 对 象 ,， 进行 中 性 检验 和 错 
配 分 布 分 析 。 结 果 显 示 为 : Тайша” s D 为 
-2.30391 ， 为 极 显 著 性 负 值 (P<0.01); Fu’s Fs 为 
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Fig. 1 








1 31 个 桔 小 实 蝇 地 理 种 群 单 倍 型 网 络 图 


Haloptype networks of 31 populations of Васігосега dorsalis 




















圆 面积 与 单 倍 型 频率 成 正比 。Circular area is proportional to the haloptype frequency. 


DCYD Сізде 1 


DCJW Clade 2 


图 2 31 个 桔 小 实 蝇 地 理 种 群 系统 发 育 树 
Fig. 2 Phylogenetic tree of 31 populations 





of Bactrocera dorsalis 


-164.545, 达到 显著 性 水 平 (P «0.02) 。 错 配 分 
析 结 果 ( 图 4) 显 示 为 单 峰 结构 。 中 性 检验 和 错 配 分 
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3 31 个 桔 小 实 晶 地 理 种 群 错 配 分 析 
Fig. 3 Mismatch distribution analysis of 31 populations 


of Bactrocera dorsalis 


析 结 果 均 表明 桂 小 实 蝇 在 中 国 地 区 曾经 存在 大 规模 
扩张 。 
2.8 桔 小 实 蝇 短片 段 COI 序列 变异 分 析 和 遗传 多 
样 性 

192 条 COI 序列 (561 bp) 共 计算 出 多 态 性 位 
点 78 个 ， 占 该 段 ОМА 序列 的 13. 90%。 其 中 自 
裔 位 点 33 个 , 占 该 段 DNA 序列 的 5. 8896 ; 简约 信 
息 位 点 45 个 , 占 该 段 DNA 序列 的 8.02% 。 共 鉴别 
出 单 倍 型 109 个 , 单 倍 型 多 样 性 为 0. 97666,， 核 背 
酸 多 样 性 为 0. 00929， 平 均 核 音 酸 差异 数 为 
5.21259, 
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到 4 23 个 中 国 桔 小 实 晶 种 群 错 配 分 析 


Fig. 4 Mismatch distribution analysis of 23 populations 


























of Bactrocera dorsalis from China 


з 讨论 


3.1 СОГ 基因 (1 496 bp) 在 种 群 遗 传 分 化 中 的 效用 

将 本 研究 获得 的 192 条 1 496 bp СОТ 序列 均 剪 
切 为 561 bp COL JE], 对 比 二 者 间 的 多 态 性 位 点 ， 
结果 显示 : 1 496 bp СОТ 序列 的 多 态 性 位 点 , 特别 
是 简约 信息 位 点 较 561 bp СОТ 序列 多 , 因而 1 496 
bp COI 序列 能 够 为 种 群 遗传 分 化 分 析 提 供 更 多 的 
言 息 位 点 。 对 比 二 者 间 的 遗传 多 样 性 指数 , 结果 显 
示 : 1 496 bp COIL 序列 除 核 背 酸 多 样 性 外 , 单 倍 型 
数 单 倍 型 多 样 性 和 平均 核 背 酸 差 异 数 均 高 于 561 
bp COI 序列 。561 bp COL 序列 的 核 苷 酸 多 样 性 
《0.00929 ) 高 于 1 496 bp COL 序列 的 核 苷 酸 多 样 性 
(0.00663), 是 由 于 561 bp COI 序列 的 总 多 态 性 位 
点 占 该 段 DNA 序列 的 比例 (13.90% ) 大 于 1 496 bp 
COI 序列 的 总 多 态 性 位 点 占 相应 РМА 序列 的 比例 
(12.90% )„ 1496 bp COI 序列 的 单 倍 型 数 和 单 倍 
型 多 样 性 均 高 于 561 bp СОТ 序列 , 说 明 单 倍 型 的 鉴 
别 因 多 态 性 位 点 的 增加 而 更 为 细 化 ,进而 使 种 群 间 
单 倍 型 分 析 更 为 精确 。 综 合 分 析 序列 多 态 性 位 点 和 
遗传 多 态 性 , 表明 СОІ 基因 (1 496 bp) 能够 为 种 群 
分 化 研究 提供 更 多 的 遗传 信息 。 
3.2 中国 桔 小 实 晶 的 扩散 路 径 

23 个 中 国 桔 小 实 晶 种 群 的 中 性 检验 和 错 配 分 
布 分 析 均 表明 : 桔 小 实 蝇 在 我 国 曾经 存在 大 规模 扩 
Ж, 而 且 中 国 原 发 生地 区 种 群 与 新 发 生地 区 种 群 间 
存在 基因 交流 或 遗传 分 化 程度 很 弱 ( 表 3)。 男 外 ， 
单 倍 型 分 析 显 示人 台北 、 香 港澳 门 . 万 宁 、 深 圳 珠海 、 
瑞丽 、 东兴 ` 长 乐 和 贵阳 等 10 个 中 国 原 发 生地 区 种 
群 与 新 发 生地 区 种 群 间 存 在 共享 单 倍 型 , 新 发 生地 
区 的 种 群 与 这 10 个 原 发 生地 区 种 群 联系 密切 。 并 
























































且 在 种 群 系统 发 育 分 析 中 , 原 发 生地 区 种 群 分 为 两 
大 分 支 , 东南 分 支 和 西南 分 支 。 通 过 本 人 研究 的 结果 
分 析 , 推断 我 国 的 桔 小 实 晶 以 东南 地 区 和 西南 地 区 
为 源头 向 内 陆 扩散 , 本 研究 结果 与 李 伟 丰 等 (2007 ) 
和 万 宣 伍 (2012 ) 的 研究 结果 一 致 。 

云南 推测 为 中 国 桔 小 实 蝎 的 最 早 人 侵 地 ( 李 伟 
丰 等 , 2007; 李云龙 等 , 2009; 吴 仲 真 等 , 2011; 万 
宣 伍 , 2012) 。 本 研究 中 , 云南 种 群 (瑞丽 和 动 脂 ) 
的 遗传 多 样 性 低 于 老挝 和 泰国 种 群 ( 沙 缴 府 和 四 色 
菊 府 ) 而 且 云 南 种 群 与 老挝 、 泰 国 种 群 间 存在 基因 
交流 或 遗传 分 化 较 弱 , 此 外 , 从 老挝 、 泰 国 种 群 到 
云南 种 群 的 迁移 率 高 于 从 云南 种 群 到 老挝 、 泰 国 种 
群 的 迁移 率 。 在 种 群 系统 发 育 树 中 (图 2), 动 脂 、 
瑞丽 与 老挝 和 泰国 种 群 属 同一 个 分 支 。 本 研究 检测 
从 东兴 珠海 种 群 到 云南 种 群 的 迁移 率 高 于 云南 种 
群 到 东兴 、 珠 海 种 群 的 迁移 率 , 表明 云南 可 能 受 老 
挝 、 泰 国 等 东南 亚 地 区 和 中 国 沿海 地 区 的 入 侵 。 男 
有 研究 结果 表明 , 云南 受到 缅甸 老挝、 泰国 、 柬 埔 
寨 等 东南 亚 地 区 和 中 国 南 海地 区 的 多 源 入 侵 (Shi ей 
al., 2012) , 本 研究 结果 与 其 一 致 。 另 外 , 云南 种 群 
与 中 国 桔 小 实 晶 种 群 间 存 在 程度 不 等 的 遗传 分 化 
(483). ПІН, 除 共享 单 倍 型 H104( 瑞丽 和 樟 木 ) ， 
云南 种 群 的 单 倍 型 均 为 种 群 特 有 单 倍 型 。 吴 仲 真 等 
(2011) 研究 中 国 及 东南 亚 地 区 30 个 桔 小 实 晶 种 群 
的 遗传 分 化 时 ,虽然 提出 云南 是 中 国 最 早 的 人 侵 
地 , 但 是 同时 也 推测 广西 为 又 一 较 早 的 人 侵 地 。 因 
此 , 我 们 推测 ,云南 不 是 中 国 桔 小 实 蝎 的 近期 
来 源 。 

珠海 东兴 长乐 和 贵阳 种 群 与 较 多 的 新 发 生地 
区 种 群 存在 基因 交流 。 在 种 群 系统 发 育 树 中 ,贵阳 
和 东兴 种 群 属于 西南 分 支 , 珠海 和 长 乐 种 群 属于 东 
南 分 支 。 此 外 ,从 珠海 种 群 到 长 乐 、 东 兴 种 群 的 迁 
移 率 高 于 从 长 乐 、 东 兴 种 群 到 珠海 的 迁移 率 ， 从 珠 
海 \ 长 乐 种 群 到 贵阳 种 群 的 迁移 率 高 于 贵阳 种 群 到 
珠海 .长 乐 种 群 的 迁移 率 , 由 此 推测 ,， 桔 小 实 蝎 和 人 
侵 到 珠海 .长 乐 等 东南 地 区 , 向 东兴 、 贵 阳 等 西南 地 
区 扩散 , 并 以 东南 地 区 和 西南 地 区 为 源头 向 中 国内 
ШИ, 其 中 珠海 长乐. 东兴 和 贵阳 种 群 为 中 国内 
陆 种 群 的 较为 有 影响 力 的 源头 。 桔 小 实 晶 新 发 生 区 
来 源 结 果 , 对 水 果 和 蔬菜 的 调运 和 检疫 进而 防 控 桔 
小 实 晶 的 传播 .扩散 有 重要 的 意义 。 
33 新 发 生地 区 桔 小 实 蝇 的 入 侵 来 源 

中 国 新 发 生地 区 种 群 与 原 发 生地 区 种 群 间 存 在 
基因 交流 或 遗传 分 化 程度 较 弱 ( 表 3) 。 另 外 , 在 种 
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群 系统 发 育 树 中 (图 2), 桔 小 实 蝇 新 发 生地 区 种 群 
在 东南 、 西 南 两 大 分 文中 均 有 分 布 。 因 此 , 新 发 生 
地 区 种 群 来 源 于 东南 地 区 和 西南 地 区 种 群 。 但 是 ， 
合肥 、 亚 东 等 种 群 的 遗传 多 样 性 高 于 东南 亚 种 群 和 
东南 沿海 种 群 ( 表 2), 因此 , 桔 小 实 蝇 新 发 生 区 可 
能 遇 到 多 个 种 群 的 多 次 人 侵 ( 万 宣 伍 , 2012) 。 本 研 
究 依 据 种 群 遗传 分 化 研究 结果 ， 可 推测 新 发 生地 区 
的 桔 小 实 晶 来 源 。 以 合肥 种 群 为 例 : 该 种 群 核 昔 酸 
多 样 性 较 云南 广东 以 及 台湾 等 原 发 生地 区 高 ( 表 
3) ,表明 合肥 种 群 可 能 有 多 个 来 源 。 合 肥 种 群 单 倍 
型 包括 H14, H22 和 H24 等 单 倍 型 ， 其 中 Н14( Ж 
X KR 香港、 衡阳 、 苏州、 赣州 上海、 深圳 和 武 
汉 ) 和 H24( 苏 州 和 上 海 ) 两 个 单 倍 型 均 被 苏州 和 上 
海 共 享 。 但 是 , 合肥 种 群 与 上 海 种 群 的 遗传 分 化 系 
数 (0. 00977) 小 于 与 苏州 种 群 的 遗传 分 化 系数 
(0.02859), 因此, 推测 安徽 合肥 种 群 的 来 源 之 一 
是 上 海 。 单 倍 型 网 络 中 , 单 倍 型 H22 和 НІ4 分 别 
衍生 于 HI17( 长 乐 和 人 台北) 和 H107( 珠 海 ) 。 合 肥 种 
群 与 长 乐 种 群 的 遗传 分 化 系数 ( - 0. 13051) 较 与 台 
北 种 群 的 遗传 分 化 系数 ( — 0. 09000) ]\, 因而 单 倍 
型 H22 衍生 于 长 乐 的 机 率 较 大 。 可 以 推测 , 合肥 村 
小 实 晶 来 源 于 长 乐 . 珠 海 和 上 海 。 同 样 ， 可 推测 樟 
木 桔 小 实 蝇 种 群 主要 来 源 于 武 隆 .贵阳 和 瑞丽 ; 武 
汉 桂 小 实 晶 种 群 主要 来 源 于 长 乐 、 珠 海 和 合肥 。 昆 
明 、 上 海 和 重庆 地 区 的 桔 小 实 晶 种 群 均 来 源 于 多 个 
种 群 (Shi et al., 2005; 万 宣 伍 等 , 2010; 吴广 超 等 ， 
2012), 与 本 文 合肥 种 群 . 樟 木 种 群 和 武汉 种 群 具有 
多 个 来 源 的 推测 结果 相似 。 

本 研究 以 CO 长 序列 对 村 小 实 蝇 种 群 进行 遗 
传 分 化 研究 , 发 现 广 东 、 福 建 广 西 和 贵州 种 群 为 中 
内 陆 种 群 较为 有 影响 力 的 源头 ; 并 根据 遗传 多 样 
TE .遗传 分 化 . 单 倍 型 分 析 , 推测 新 发 生地 区 桔 小 实 
晶 的 来 源 , 例如 合肥 种 群 主要 来 源 于 长 乐 .珠海 和 
上 海 , 但 本 研究 结果 表明 村 小 实 蝇 地 理 种 群 间 的 遗 
传 分 化 程度 较 弱 。 在 遗传 分 化 程度 弱 时 , 需要 更 多 
的 样本 量 或 基因 来 保证 分 析 的 可 靠 性 ( Kalinowski， 
2005) 。 因 此 , 下 一 步 研究 还 需 增加 不 同 地 理 种 群 
的 样本 量 或 基因 数目 , 桔 小 实 晶 种 群 遗 传 结构 分 析 
结果 才 会 更 可 靠 。 
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